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Abstract We study pseudogap behaviors of ultracold Fermi gases in the BCS-BEC 
crossover region. We calculate the density of states (DOS), as well as the single-particle 
spectral weight, above the superfluid transition temperature Tc including pairing fluc- 
tuations within a T-matrix approximation. We find that DOS exhibits a pseudogap 
structure in the BCS-BEC crossover region, which is most remarkable near the uni- 
tarity limit. We determine the pseudogap temperature T* at which the pseudogap 
structure in DOS disappears. We also introduce another temperature T** at which 
the BCS-like double-peak structure disappears in the spectral weight. While one finds 
T* > T** in the BCS regime, T** becomes higher than T* in the crossover and BEC 
regime. We also determine the pseudogap region in the phase diagram in terms of 
temperature and pairing interaction. 
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1 Introduction 

Recently, the BCS-BEC crossover has been realized in ultracold Fermi gases[T]. In this 
phenomenon, using a tunable pairing interaction associated with a Feshbach resonance, 
one can study Fermi superfiuids from the weak-coupling BCS regime to the strong 
coupling BEC regime in a uinfied manner. Because of this advantage, superfiuid Fermi 
gases would be also useful for the study of high-Tc cuprates with a strong pairing 
interaction. 

In the under-doped regime of high-Tc cuprates, the so-called pseudogap structure 
has been observed in the density of states (DOS) [2]. As the origin of the pseudogap. 
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strong pairing fluctuations has been proposed[3]|4l[5l[6]. However, because of tlie com- 
plexity of this system due to strongly correlated electrons, other possibilities, such as 
antiferromagnetic spin fluctuations and a hidden ordered state, have been also dis- 
cussed. Thus, to confirm the pairing fluctuation scenario, another simple system only 
having superfluid fluctuations would be useful. 

The cold Fermi gas system meets this demand. It is much simpler than high- 
Tc cuprates, and pairing fluctuations dominate over the BCS-BEC crossover physics. 
Indeed, the pseudogap phenomenon in this system has been recently predicted [3l[6l 
[7]. Although the s-wave pairing symmetry of superfluid Fermi gas is different from 
the d-wave one in high-Tc cuprates, we can still expect that the study of pseudogap 
phenomenon in cold Fermi gases would be helpful in understanding the under-doped 
regime of high-Tc cuprates. Since a photoemission-type experiment has recently become 
possible in cold Fermi gases [8], observation of strong-coupling effects on single-particle 
excitations is now possible within the current technology. 

In this paper, we investigate pseudogap behaviors of atomic Fermi gases above the 
superfluid transition temperature Tc- Including pairing fluctuations within a T-matrix 
approximation, we calculate DOS and single-particle spectral weight, we examine how 
pairing fluctuations affect them over the entire BCS-BEC crossover region. We also 
determine the pseudogap regime in the phase diagram in terms of temperature and the 
strength of pairing interaction. 



2 Formalism 

We consider a uniform two-component Fermi gas described by pseudospin a —^,1- 
For a broad Feshbach resonance (which all the current experiments are using), it is 
known that one can safely study the interesting BCS-BEC crossover physics by using 
the ordinary BCS model[T], given by 

H ^^^pcl'rCpa - t^X!X!'^P+<j/2T^-p+q/2i^-P'+'3/2iV+q/2T- W 
p,cr q p^p' 

Here, Cpa is an annihilation operator of a Fermi atom with pseudospin u |. = 
Ep — ^ — /2m ~ II is the kinetic energy, measured from the chemical potential ji, 
where m is an atomic mass. The pairing interaction {U > 0) is assumed to be tunable 
by a Feshbach resonance. In cold atom physics, the strength of pairing interaction is 
conveniently described in terms of the parameter {kpas)~^ , where as is the s-wave 
scattering length and kp the Fermi momentum. In this scale, the BCS limit and BEC 
limit are, respectively, given by (fepas)^^ <Si —1 and (fcpas)^^ ^ +1. The region 
— 1 <i (fcptts)"^ ^ -I- 1 is referred to as the crossover region. The relation between U 
and as is given by [9] A-Kas/m = —U/[l — U '^^{m/p'^)]. 

The single-particle thermal Green's function is given by Gp{iLOn) = 1 / {{Gp{iujn))^^ — 
E{p,iujn))- Here, ujn is the fermion Matsubara frequency, and Gp{iLJn) ~ l/{iu!n —^p) 
is the non-interacting Fermi Green's function. The self-energy S{p, iuJn) involves effects 
of pairing fluctuations. In this paper, we include strong-coupling corrections within the 
T-matrix approximation [3]|4l[6]. The resulting self-energy has the form 

E{p,iLOn) = T ^ T(q,i!/„)Gq-p(ii^n - «'^n)e**''^"~""''*i (2) 
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where Vn is the boson Matsubara frequency. The particle-particle scattering matrix F 
(describing pairing fluctuations) is given by r{q,ivn) ~ —U/[l — Un[q,ivn)], where 
n{q,ivn) = Tj2p.^,^ G'p+q/2(*'^» + «'^»)G"p+q/2(~*^») a pair-propagator. 

To discuss the pseudogap phenomenon above Tc, we need to determine Tq. Following 
Ref.[TO], we employ the Thouless criterion [TT]. r{q = 0,iun ^ 0,T = Tc)"^ = 0, and 
solve this equation, together with the equation for the number of fermions, 

iV = 2r ^ e"^"^Gp{iiJn). (3) 

The above treatment can describe the smooth crossover behavior of Tc and ^ in the 
BCS-BEC crossover [6]|10|. Namely, starting from the weak-coupling BCS regime, Tc 
gradually deviates from the mean-field result to approach Tc — 0.218eF of an N/2 ideal 
molecular Bose gas (where ep is the Fermi energy). The chemical potential monoton- 
ically decreases from £p in the crossover regime to be negative in the BEC regime 
{{kpas)~^ > 0.35). The negative /i indicates the formation of two-body bound states, 
so that the BEC regime is well described by a molecular Bose gas, as expected. 

Above Tc, we solve Eq. @ to determine fi. DOS p{uj) and the spectral weight 
A{p, u) are, respectively, evaluated from the analytic continued Green's function, as 

p{oj) = lmGp{iLOn UJ+), (4) 

V 

A{p,u) = -—Gp{iujn (5) 
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Fig. 2 Temperature dependence of the density of states p{i^). To clearly show how the pseu- 
dogap disappears, we have offset results for T > Tc. The horizontal line in the left end of each 
curve indicates zero density of states. 



3 Results 

Figure [T] shows DOS at Tc. In the BCS side ((fcpas)""^ < 0) (panel (a)), we find a 
pseudogap (dip) structure around u — 0, which evolves as one approaches the unitarity 
limit {{kpas)~^ = 0). Since we only include pairing fluctuations, this pseudogap purely 
originates from pre-formed pairs in the normal state. 

In the BEC regime (where fi < 0), when we only include the negative fi and 
ignore other strong coupling effects, DOS has a finite gap as p{lo) oc lu + Indeed, 
Fig[ljb) shows that DOS continuously changes into the fully gapped structure. Since 
2\p\ equals the binding energy of a molecule in the BEC limit 10 , the almost fully 
gapped structure at (fcpas)^^ = 1 in Fig[ljb) indicates that the system is rather close 
to a molecular Bose gas than a Fermi atom gas. 

Figure [2] shows DOS above Tc. As expected, the pseudogap gradually disappears as 
one increases the temperature. When we define the pseudogap temperature T* at which 
the dip structure in DOS disappears, we obtain the phase diagram Fig|3] Although Tc 
calculated in the mean-field theory is usually considered as a characteristic temperature 
at which preformed pairs appear, Fig|3] shows that the pseudogap temperature T* 
evaluated from DOS is actually much lower than the mean-field value of Tc {Tbcs in 
this figure). 

In the BCS theory, the spectral weight below Tc has the double-peak structure as 



A{p,Lj) 



5{uj-Ep) + - 



En 



S{uj + Bp) 



(6) 
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Fig. 3 Pseudogap temperatures T* and T** determined from DOS and spectral weight, re- 
spectively. In this figure, SF is the superfluid phase, PG is the pseudogap regime, and MB is 
the molecular Bose gas regime. Above T* or T**, the system is regarded as a normal Fermi 
gas (NF) with weak pairing fluctuations. 



Here, Ep = -y/^p + is the Bogoliubov excitation spectrum, where A is the superfluid 
order parameter. The minimum value of peak-to-peak energy equals 2A. In the simple 
mean-field theory, this double-peak structure only exists below Tc- However, when one 
includes pairing fluctuations, the double-peak structure still remains above Tc[3]l6], as 
shown in Fig|4] This characteristic structure disappears at a certain temperature T** 
(See Fig|4ja) and (b)), which we can define as another pseudogap temperature. 

When we plot T** in FigE] we find that it does not coincide with T* . While the 
latter is higher in the BCS side, one finds T** > T* in the EEC side. In the BCS 
side, although the double-peak structure is absent when T > T** , pairing fluctuations 
still strongly affect A{p,uj) around the Fermi level p ~ ^/2mfi, leading to a broad and 
low single peak structure. Then, since DOS is given by the momentum summation of 
A{p,Lj), DOS around a; ~ (which is dominated by A{p ~ ^2m/7, oj)) is suppressed, 
giving the dip structure. On the other hand, as shown in Fig|4fc), the lower peak in 
A{p,u!) soon becomes broad above Tc in the BEG regime, so that the effect of lower 
peak is easily smeared out in the momentum summation in Eq. ((SJ. As a result, DOS 
does not reflect the double-peak structure in A{p,u!) when T* < T < T** . 

The above result indicates that the pseudogap temperature depends on what we 
measure. When we measure a quantity where DOS is crucial, such as the specific heat, 
T* would be the crossover temperature between the pseudogap regime (PG) and normal 
Fermi gas regime (NF). On the other hand, when we consider a quantity dominated by 
the spectral weight, such as the photoemission-type experiment done by JILA group [S], 
T** would work as the crossover temperature between PG and NF. 

As discussed previously, the pseudogapped DOS continuously changes into fully 
gapped one in the BEG regime, reflecting that the system reduces to an N/2 molecular 
Bose gas (MB). As a crossover temperature between PG and MB regime, the molecular 
binding energy Ei^i^^j^ would be useful, because thermal dissociation of molecules is 
suppressed when T <, -Ebind- In the present case, noting that iSbind — 1 2^1 in the BEG 
regime (where jj, < 0), one may conveniently determine the pseudogap regime as the 
region surrounded by Tc, 2\iJ,\, and T* or T** in Fig|3] 
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Fig. 4 Spectral weight A{p,u!) as a function of uj. In each panel, we take the momentum 
where the peak-to-peak energy is minimum: (a) p/kp = 0.91, (b) 0.83, and (c) 0.01. 



4 Summary 

To summarize, we have discussed effects of pairing fluctuations on single-particle prop- 
erties of a cold Fermi gas above Tc. Within the framework of T- matrix approximation, 
we showed how the pseudogap appears in the density of states and the single-particle 
spectral weight over the entire BCS-BEC crossover region. We have also determined 
the pseudogap regime in the BCS-BEC crossover phase diagram. This phase diagram 
would be useful for the observation of pseudogap phenomenon in this system. 
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